Impurity-induced bound states in superconductors with topological order 



Fei Wang Qin Liu Tianxing Ma '•^•*, and Xunya Jiang ^ 
' Department of Physics, Beijing Normal University, Beijing 100875, China 
^ State Key Laboratory of Functional Materials for Informatics, 
Shanghai Institute of Microsystem and Information Technology, CAS, Shanghai 200050, China and 
' Beijing Computational Science Research Center, Beijing 100084, China 
(Dated: October 12, 2012) 

The study of classical spins in topological insulators [Phys. Rev. B 80, 115216 (2009)] is generalized to 
topological superconductors. Based on the characteristic features of the so-called F-function, Bogoliubov- 
de Gennes Hamiltonian for superconductors is classified to positive, negative, and zero "gap" categories for 
topologically trivial and nontrivial phases of a topological superconductor as well as a BCS superconductor 
respectively. It is found that the F-function determines directly the presence or absence of localized excited 
states, induced by bulk classical spins and nonmagnetic impurities, in superconducting gap and their persistence 
with respect to impurity strength. Our results provide an alternative way to identify topologically insulating and 
superconducting phases in experiments while without resorting to the surface properties. 
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I. INTRODUCTION 

Time-reversal (TR) invariant topological insulators (TI) are 
new state of matter in condensed matter physics, which have 
a full insulating gap in the bulk, but gapless edge or surface 
states consisting of an odd number of Dirac fermions. ' Materi- 
als for TR invariant TIs, including HgTe/CdTe quantum wells, 
Bii_vSbt alloy, binary (Bi2Te3, Sb2Te3, Bi2Se3) and ternary 
(TlBiTe2, TlBiSe2) compounds, have been confirmed exper- 
imentally according to the above characteristic definition, by 
either transport measurements of the quantized edge conduc- 
tance 2e^/li, or angle-resolved photoemission spectroscopy 
measurements of the linear dispersion of the surface states as 
well as its odd numbers^i However, all these measurements 
aim only at the surface properties of the TR invariant TIs. 
While in a recent worfc2. by two of the authors, through ex- 
amining the localized excited states (LES) in the insulating 
gap induced by bulk impurities in quantum anomalous Hall 
(QAH) systems, in two-dimensional (2D) HgTe/CdTe quan- 
tum spin Hall (QSH) systems^i as well as in TR invariant 3D 
strong TIs^^ they point out that distinctive behaviors between 
topological and conventional insulators exist, which help to 
distinguish topologically different insulating phases experi- 
mentally without resorting to surface features. 

Soon after the discovery of TR invariant TIs, the study was 
generalized to TR invariant and breaking topological super- 
conductors (TSC)ji'^ This generalization is natural because 
there is a direct analogy between superconductors and insu- 
lators where the Bogoliubov-de Gennes (BdG) Hamiltonian 
for the quasiparticles of a superconductor is analogous to the 
Hamiltonian of a band insulator, with the superconducting gap 
corresponding to the band gap of the insulator. Instead of 
Dirac fermions, the gapless surface states of a TR invariant 
TSC consist of odd number of Majorana cones, which only 
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have half the degree of freedom of Dirac fermions. While 
for TR breaking TSC, the vortex core of which carries an 
odd number of Majorana zero modes,' giving rise to non- 
Abelian statistics and providing a possible platform for topo- 
logical quantum computing^ii Several ways to realize topo- 
logical superconductivity by making use of the superconduct- 
ing proximity effect on the 2D surface states of 3D TR invari- 
ant TIs,— or on the 2D TR breaking TIs,— and on semicon- 
ductors with strong Rashba spin-orbit couplingi^ have been 
imposed. However, no definitive proofs in experiments have 
yet been found so far. 

Gaped systems can be classified into ten symmetry classes, 
among which four are the BdG classes for superconductorsi^ 
Bearing in mind the extreme analogy between the BdG Hamil- 
tonian of a superconductor and that of a band insulator, we 
show in this work that the BdG Hamiltonian for supercon- 
ductors can also be assorted in another viewpoint into three 
categories of positive, negative, and zero "gap" superconduc- 
tors according to the so-called F-functionsi2iii By doping into 
bulk magnetic (which we treat as classical spinsiS) and non- 
magnetic impurities, we show that the characteristic proper- 
ties of the F-functions determine directly the presence or ab- 
sence of LES, induced by the impurities, in superconducting 
or insulating gap, as well as the persistence of the LES with 
respect to impurity strength. Based on this observation, a 
generic method is proposed, by testing the response to clas- 
sical spins as well as nonmagnetic impurities in the bulk, to 
differentiate the topologically nontrivial superconducting or 
insulating phases from the trivial ones in general, and also to 
distinguish a TSC from a conventional BCS superconductor 
in particular 

Specifically speaking, using the T-matrix methodr^ it is 
found that, similarly as in TIs? in TSCs there are always four 
LES in the superconducting gap for spin-dependent potential, 
whereas only two such LES for ordinary potential. Moreover, 
these LES survive under arbitrary impurity strength. While 
when the TSC transits into the topologically trivial phase, the 
LES exist only at small impurity strength but then disappear 
into bulk bands at strong impurity strength. The classical spins 
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FIG. 1; (Color online) Energy spectra of rf + id superconductors 
with open boundary condition in x-direction. Parameters are taken 
as Aj2_,,2 = 2, A„ = -2, and all energies are in unit of t. (a) Spec- 
trum for a topologically trivial phase at yu = -3. (b) Spectrum for a 
topologically nontrivial phase at yu = 3. Two gapless edge modes in 
the superconducting gap are marked in red. 



II. BRIEF REVIEW OF SYMMETRY CLASSIFICATION 
OF BDG HAMILTONIAN FOR SUPERCONDUCTORS 

In this section, we briefly review the symmetry classifica- 
tion of BdG Hamiltonian of superconductors. It is concluded 
that within all topologically nontrivial categories of the four 
BdG classes in dimension two and three,— only class C in di- 
mension two is left unaddressed for a complete discussion. 

In a recent pioneering work by Schnyder et air- ten sym- 
metry classes of single-particle Hamiltonians for gapped sys- 
tems are classified according to the presence or absence of 
TR symmetry, PHS, as well as sublattice symmetry. Among 
which, four symmetry classes (named as D, Dili, C, CI) arise 
in the BdG Hamiltonian for superconductors because of the 
definitive PHS but the alternative presence of TR symmetry 
or SU(2) spin rotational symmetry. 

A general form of BdG Hamiltonian for the dynamics of 
quasiparticles deep inside the superconducting state of a su- 
perconductor can be written asi21 











-t) 


CD 



(1) 



localized in a BCS superconductor have been discussed by H. 
Shiba in 196^ In the viewpoint of F-functions, BCS super- 
conductors fall into the critical category between the topolog- 
ically trivial and nontrivial phases of a TSC. It is shown that 
two LES appear only for spin-dependent potential but no such 
LES exist for ordinary potentials. Moreover, these two LES 
persistently tend nearer to the band edges at strong impurity 
strength, which is in a tricky contrast to the trivial phase of a 
TSC. Therefore, through the observation of the LES in super- 
conducting (band) gap induced by bulk impurities, the poten- 
tial strength of both spin-dependent and ordinary impurities 
can be used to tell the topologically nontrivial from the trivial 
phases of a TSC (TI), whereas nonmagnetic impurities can be 
used to distinguish a TSC from an ordinary (BCS) supercon- 
ductor. 

The rest of this paper is organized as the following. In Sec- 
tion [n] we first briefly review the symmetry classification of 
BdG Hamiltonian for superconductors.^ Through this classi- 
fication, all categories of BdG Hamiltonian can be reduced, 
if the additional particle-hole symmetry (PHS) is imposed, 
to the model Hamiltonian of either the 2D QSH system^i or 
the 3D strong TI,^ or half the spin degree of freedom of the 
above two which breaks the TR symmetry i^i^S The responses 
to classical spins and nonmagnetic impurities of all these cor- 
respondences in TIs have been studied in Refi^, with the only 
exception being class C superconductors In Section|III] we 
study the classical spins in ^^^2_,,2 + idxy superconductors as 
an specific example of class C. It is shown that its response 
to both spin-dependent and ordinary potentials obey the same 
rules as those^ in TIs. Therefore, observing the LES in super- 
conducting (band) gap induced by bulk impurities serves as 
an effective criterion to tell the topological superconducting 
(insulating) phase. This work is finally concluded in Section 
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where Ck and c| can be either column or row vectors, eu - 
12m - i^is the single-particle energy dispersion, Ak - (dk ■ 
£r)(/cr,,) is the gap parameter with cr being the electron spin 
Pauli matrix vector, and is a 3D vector in spin space as a 
function of momentum k. 

In class D, where neither SU(2) invariance nor TR symme- 
try is present, there exist topologically nontrivial supercon- 
ducting phases in both ID and 2D. A typical example of BdG 
Hamiltonian in this class in 2D is the spinless chiral p-wave 
ip+ip) superconductors, where the gap parameter is explicitly 
A*; - K{k^ - iky), A e M. We notice that this BdG Hamilto- 
nian of TR breaking p-wave superconductors is nothing else 
but the BHZ model^ of a QSH system with half the spin de- 
gree of freedom, or the QAH system studied by Qi et air- 
by dropping the terms proportional to the identity matrix to 
maintain the generic PHS. In this sense, a TSC can be viewed 
as a TI with PHS. Classical spins in QSH and QAH systems 
have been studied in the context of TI,i where the term which 
breaks the PHS is shown to be unimportant to the topological 
properties of the systems under study. Therefore it is safe to 
conclude that the p + ip superconductors have nontrivial re- 
sponse to spin-dependent (ordinary) potential, where there are 
always two (one) LES in the superconducting gap in topolog- 
ically nontrivial phase, which disappear in trivial phase. 

In class Dili, though full SU(2) invariance is still absent but 
TR symmetry is restored. The simplest way to regain TR sym- 
metry is to make two copies of the TR breaking system but 
with the whole being a TR conjugate pair. In this class, topo- 
logically nontrivial superconducting phases exist in all three 
dimensions. An interesting example in 2D is the equal super- 
position of two chiral p-wave superconductors with opposite 
chiralities (p^ + ipy and p^ - //?,, waves), namely the helical 
p + ip superconductors, where the d-vector is explicitly dk - 
A{-k„ ky, 0). In the basis of (c^ c_*) = (c^^ c-ki c|j^ c_y), we 
see that the BdG Hamiltonian of the helical p + ip supercon- 
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ductors is completely identical to that of the BHZ model in 
HgTe/CdTe quantum wells with PHS, where spin up (down) 
electrons form + ipy (p^ - ip^.) Cooper pairs. Hence the 
response of such helical p + ip superconductors is similar to 
that of the QSH systems, where persistent LES exist in the 
superconducting gap. In 3D, a member of class Dili is the 
Balian-Werthamer (BW) stated of the B phase of liquid ^^He 
described by the d-vector dt = A(^j, ky, k.). In the same basis 
as the helical p + ip superconductors, the model Hamiltonian 
of the BW state is the same as that of a 3D TI, which reduces 
exactly to the 2D QSH system in the limit k- ~ (k^) ^ 0. 
We have also studied the sp-d exchange coupling in 3D TIs^^ 
again the LES behave distinctively in topologically nontriv- 
ial and trivial phases, which is also valid to the BW state of 
topological superfluids. 

With the presence of full SU(2) invariance, the BdG Hamil- 
tonians fall into C and CI classes in the absence and presence 
of TR symmetry. Nontrivial topological phases exist in 2D in 
class C, while in 3D in class CI. A full focus on the five cate- 
gories of TR invariant TIs in 3D out of ten symmetry classes 
have been laid on in the pioneering work by Schnyder et al.r- 
therefore for a complete discussion the only unaddressed one 
is the 2D case in class C, which is our main attention in this 
work. The interesting example in class C in 2D is the TR 
breaking superconductors with (d + ;(i)-pairingiai Under the 
basis (c^ C-k) - (cl^^ c-ki), the gap parameter ind + id super- 
conductors is 



At = Aj,2 ^yiik;. -ky) + ihxyk^ky 



(2) 



where A j-2_,,2 and An are real amplitudes. 

In the following, we study the topological properties of 
the d + id superconductors as well as its response to clas- 
sical spins and nonmagnetic impurities in bulk using the T- 
matrix method. It is found that the d + id superconductors 
show the same behavior as that of TIs. Specifically, in topo- 
logically nontrivial phase, there are four LES in the super- 
conducting gap for magnetic impurities and two LES for non- 
magnetic impurities at arbitrary impurity strength. These LES 
then disappear into the bulk bands at large impurity strength 
as the system transits into the topologically trivial phase. For 
comparison, no LES appear with nonmagnetic impurities in 
a BCS superconductor Therefore, by testing the response to 
both magnetic and nonmagnetic impurities, not only a topo- 
logically nontrivial phase can be distinguished from a trivial 
one, but also a conventional BCS superconductor can be dis- 
criminated from a TSC. 



III. CLASSICAL SPINS IN d^i 
SUPERCONDUCTORS 
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In this section, we study d + id superconductors described 
by Eq.© as an specific example in class C in parallel to the 
previous work in Ref.\ The results are stated in the language 
of superconductors in particular, which are also valid for in- 
sulating systems in general. 

The single-particle Hamiltonian of d + id superconductors 
in momentum space is written as h{V) - ^,(k)cr", where a - 
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FIG. 2: (Color online) F-functions versus w which lies in bulk super- 
conducting gap for different /j's. (a)-(d) F-functions in topologically 
nontrivial phase, (e)-(h) F-functions in topologically trivial phase. 
The electron band F-function Fg(oS) is shown in red lines while the 
hole band F-function Fi,(a>) is shown as green lines. Parameters are 
taken as sgn(A,.2_,,2 ■ A„.) = -1 and all energies are measured in units 
off. 



1,2,3, and in tight-binding model, f(k) = (-Aj.2_,,2(cos - 
cos ky), -A„. sin k^ sin ky, f(4 - cos k^ - cos ky) - jj), with t 
being the hopping energy. The Hall conductivity of this two- 
band system at zero temperature when the chemical potential 
lies inside the band gap is calculated using the standard Kubo 
formula as^ 



dk^dkyi ■ dki x dkX 



(3) 



where /»(k) - fa(k)/f{k) is the unit vector along the direction 
/„ and /(k) = -y//„(k)/"(k). This Hall conductivity is related 
with the so-called Skyrmion numberSiSS by cr^y - -Q^i^y/ln, 
which is an integer, and ind + id superconductors it is shown 
to be 



2sky - 



-2sgn(A^2_^2/A;,j.), < A< < 8, 
0, p <0 or jj> 



(4) 



Therefore the bulk-edge correspondences tell us that there 
should localize two edge states at each boundary in the topo- 
logically nontrivial phase, which is indeed true as seen in 
FiglT]where typical energy spectra of d + id superconductors 
are shown. 
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To consider the single-particle excitations in the supercon- 
ducting gap, our starting point is to rewrite the BdG Hamilto- 
nian of d + id superconductors in Nambu spac( 
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where the gap parameter A^^ is given in Eq.©, and = 

{c\^ C-ki c-k^)- This is because the Green's function (GF) 
formulation in Nambu space can be easily generalized to in- 
clude the paramagnetic impurities as well as the Kondo effect. 
Also notice that the above BdG Hamiltonian of a d + id su- 
perconductor in Nambu space is formally similar to the 2D 
Luttinger model in the hole bands of a semiconductor.- In 
the presence of a localized spin or an ordinary potential in 
bulk, their interactions between conduction electrons are re- 
spectively 



//ex = (J/2) 2 C|S (T, //' = (V/2) 2 C; 

kk' kk' 
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where S = (Sx,Sy,S.) is the vector of a localized spin with 
modulus - SaS". Again in Nambu space they take the 
form 
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FIG. 3: (Color online) Real space diagonalization of the total Hamil- 
tonian //BdG + //cx in (a)-(b) and Hbag + H' in (c)-(f) as a function of 
impurity strength. The bulk states are represented by blue lines and 
LES in superconducting gap are denoted in red lines. Parameters are 
taken as sgn(A,.2_,2/A„) = -1 in all the figures, yu = ±3 in topo- 
logically nontrivial and trivial phases, and all energies are measured 
in units of t. (a)-(b) Topologically trivial and nontrivial phases for 
spin-dependent potential, (c)-(d) Topologically trivial and nontrivial 
phases with V > for ordinary potential, (e)-(f) Topologically trivial 
and nontrivial phases with K < for ordinary potential. 
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where S ± = S x ± iSy Then the full GF of ad + id supercon- 
ductor with a classical spin or an ordinary potential is obtained 
through the equation of motion methodic of GF by 



(9) 



where G°(w) = l/(<y - //Bdc) is the free GF of d + id su- 
perconductors, while the T-matrix for a localized spin and an 
ordinary potential are separately 



(JS/2YF(oj) + 



fex(w) = 



S+T- -5- 



1 - [( JS/2)Fica)r~ 



H' 



1 - H'F(aj) 



(10) 
(11) 



In the above, the T-matrix depends only on the energy be- 
cause the scattering in Eqs. (|7]i and (|8]l are both momentum- 
independent, the Ta's are Pauli matrix in particle-hole space, 
and we have defined the F-function as 

F(co) = I 2 Glioj) = diag(F, Fh F, F,,) (12) 

k 



with the diagonal elements 
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It is interesting to notice that this /"-function in d + id su- 
perconductors is diagonal too as that in the QSH system in 
HgTe/CdTe quantum wells with p-wave symmetry studied 
beforei^ This is because the off-diagonal terms are momen- 
tum integrations of the gap parameter A^., although Aa has 
quadratic dependence on momentum, the integration of real 
part cancels exactly since the integrand is symmetric with re- 
spect to kx and ky, while the imaginary part vanishes uniformly 
for that the integrand is an odd function of k^ and ky. There- 
fore, the eigen-energies for the LES induced by classical spins 
and ordinary potentials are obtained by finding the poles of the 
above T-matrix in (fTOl i and ( fTTT i respectively for energies in the 
superconducting gap as 

(JS/2)Feu,y{oj) = +1, (14) 
(V/2)F,(aj) = 1, and (V/2)F,,(aj) = -1, (15) 

which consist of four conditions when the impurity is spin- 
dependent while only two conditions when the impurity is 
spin-independent. 

To see how these conditions are satisfied for general impu- 
rity strength J and V, we numerically examine the momen- 
tum integrations in Eq.(fT3Tl in both topologically nontrivial 
and trivial phases which are identified in Eq.® by constrain- 
ing the frequency co in the superconducting gap only. Using 



the tight-binding model where — t{4 — cos A;^ - cos k^) — ji, 

|A|(:|- = A^2_ .2(cos A:^ - cosA,.)^ + sin^Aj. sin^A,., and re- 

X y ■ ■ y 

placing {l/N)Yjk ~* (l/4-7p-) Jf dkjcdky, several representa- 
tive F-functions are plotted versus frequency for topologically 
nontrivial phases in Figs|2a)-(d), and for topologically triv- 
ial phases in Figs|2le)-(h). A striking difference of the F- 
functions between topologically nontrivial and trivial phases 
is immediately recognized in these figures. First of all, it is 
seen that each of the electron band F-function Feico) and the 
hole band F-function FhiuS) winding from -oo to oo within 
the entire range of u in the superconducting gap, and there 
is a negative "gap" between Fe{a>) and in topologically 

nontrivial phase [See Figs|2a)-(d)]. Therefore, all of the four 
resonant conditions in Eq.(fT4li for classical spins can be satis- 
fied at any impurity strength JS (for both / > and / < 0), 
and we expect four persistent LES exist in the superconduct- 
ing gap. Similarly for ordinary potential, both of the resonant 
conditions in Eq.dTSll can also be met with arbitrary V (> or 
< 0), hence two persistent LES in the superconducting gap is 
expected when the impurity is spin-independent. While in a 
sharp contrast, in topologically trivial phase, both of the elec- 
tron band F-function and the hole band F-function terminate 
at some finite value in one end of number- axis whereas di- 
verge in the other end, which are separated by a positive "gap" 
[See Figs|2je)-(h)]. Consequently, only one condition in each 
of Fe{u) and in Eq.(fT4li becomes true at small JS but 

fails for strong enough impurity strength in the case of clas- 
sical spins. The case is tricky for ordinary potentials, which 
depends on the sign of product of V and //, sgn(y//). Since 
the electron band F-function F^ioS) is negative definite and 
the hole band F-function Fi,{a)) is positive definite for co in 
the superconducting gap when fi < [See Figs|2je)-(f)], no 
LES exist for repulsive interactions but two LES appear for at- 
tractive interactions. And in attractive situation, the two LES 
finally disappear at large V. Vice versa when ;U > 8. 

As an independent examination of the above statements an- 
alyzed solely from the F-function behaviors, we have also di- 
rectly diagonalized the total Hamiltonian H - H^dc+HnxiH') 
in real space as a function of impurity strength \JS \ (\V\) in 
Fig|3] In Figs|3la) and (b), real space diagonalization with 
spin-dependent potential is shown, where we see that in the 
topologically trivial phase (up panel), there are indeed two 
LES winding through the superconducting gap at small \JS |, 
which disappear at strong \ JS \. Nevertheless, four LES persist 
in the superconducting gap at arbitrary impurity strength in 
topologically nontrivial phase (lower panel), and the stronger 
the \ JS I is the deeper the LES are localized in the gap. While 
in Figs|3jc)-(f), real space diagonalization with ordinary po- 
tential is presented. We see that in topologically nontrivial 
phases for both sign of V [Figsl3ld) and (f)], two persis- 
tent LES exist in the superconducting gap. And in topolog- 
ically trivial phases, no LES appear in the superconducting 
gap when V/i < [Figl3jc)] but two impersistent LES ex- 
ist at small \V\ when V;U > [FigOIe)]. These pictures sup- 
port in perfect to our predictions merely from the F-function 
properties given in Fig|2] In addition, we could expand the 
5-function of the impurity potential and consider an a finite 
range spin-dependent scattering with V - VoF/ (jr - ro|^ + F^). 
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FIG. 4: (a) F-function in BCS superconductors where the electron 
and hole band F-functions terminate at the same value (see the dotted 
line) and form a zero "gap", (b) Eigen-energies of two LES induced 
by classical spins versus impurity strength in a BCS superconductor. 



For an extended impurity with small Vq, the impurity eff'ect 
is restricted on the impurity site as that in the localized im- 
purity case^i. For a relative large Vq, the impurity scattering 
is strongest at the impurity site and decay quickly as the lat- 
tice spacings increase, which is not important beyond few lat- 
tice spacings from the impurity citeM. Then, it is expected 
that similar results could be addressed for a finite range spin- 
dependent scattering, and detail results should be present in 
our future work. 

For comparison, we now discuss the F-functions as well 
as LES induced by classical spins and nonmagnetic impu- 
rities in BCS superconductors, which have been studied in 
early times by H. Shiba.J^ In a BCS superconductor, the cor- 
responding F-function is analytically obtained as Fbcs(w) = 
-.TrA^f ((w/Ao) + 1)/ Vl - (w/Ao)2, where Aq is the real gap 
parameter and Np is the density of state at Fermi energy in 
normal state. This F-function Fg^daij/inNf) versus energy 
is plotted in FigHJa) for |w/Ao| < 1, where we see that both 
of the electron and hole band F-functions diverge at one end 
while terminate at the same value at the other end. There- 
fore, in this sense, a BCS superconductor can be viewed as 
a zero "gap" superconductor compared to the positive "gap" 
topologically trivial and negative "gap" topologically nontriv- 
ial phases of a TSC. For classical spins in a BCS superconduc- 
tor, there are always two LES in the superconducting gap at 
energies +cjbcs = +Ao(l -(/^TrA^f /2)2)/(l -H(y57rA^f /2)2) as 
shown in FigHJb), where we see that they behave critically be- 
tween the topologically trivial and nontrivial phases exhibited 
in Fig|3] The two LES in BCS superconducting gap persis- 
tently tend nearer to the band edges as the impurity strength 
goes to infinity but never disappear. While in contrast, the two 
LES of the trivial phase in a TSC goes to the band edges as 
the increase of impurity strength and then disappear into the 
bulk at some finite value of |75 j or jVj. For ordinary potential 
in BCS superconductors, it is shown that LES never appear 
because the corresponding T-matrix has no singularities in the 
superconducting gap. 

The response of a TSC to classical spins and nonmagnetic 
impurities compared with that of a conventional BCS super- 
conductor is summarized in Table. U where those persistent 
LES which survive at arbitrary impurity strength are num- 
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bered in bold. This result provides an effective way to distin- 
guish the topologically trivial and nontrivial phases of a TSC, 
and also to tell a conventional BCS superconductor from a po- 
tential TSC, while without resorting to the surface properties. 
This is the main result of our work. 

TABLE I: The number of LES in the superconducting gap induced 
by classical spins and nonmagnetic impurities in a TSC as well as in 
a conventional BCS superconductor. The bold numbers indicate that 
the corresponding LES exist at arbitrary impurity strength. 





trivial-TSC 


nontrivial-TSC 


BCS 


75/2 


2 


4 


2 


K/2 


0, sgn(Vyu) < 


2 







2, sgn(Vyu) > 







IV. CONCLUSIONS AND DISCUSSIONS 

In conclusion, classical spins in generic TIs and TSCs are 
studied in an earlier^ and the present work. In particular, an F- 
function is defined in such systems, where the differences of 
which from electron and hole band classify BdG Hamiltonian 
of superconductors into positive, negative, and zero "gap" 
categories respectively for topologically trivial and nontrivial 



phases of a TSC and a BCS superconductor. The character- 
istic features of F-functions determine directly the presence 
or absence of LES, induced by bulk classical spins and non- 
magnetic impurities, in superconducting gap as well as their 
persistence with respect to impurity strength. The responses 
of TSC and BCS superconductors to bulk classical spins and 
nonmagnetic impurities are summarized, where it is shown 
that the potential strength of both spin-dependent and ordi- 
nary impurities can be used to tell the topologically nontrivial 
from the trivial phases of a TSC, whereas nonmagnetic impu- 
rities can be used to distinguish a TSC from an BCS super- 
conductor Our results provide an alternative way to identify 
topologically insulating and superconducting phases in exper- 
iments without resorting to the surface properties. 
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